Density Fluctuations in Molten Lithium: 
Inelastic X-Ray Scattering Study. 
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New inelastic X-ray scattering experiments have been performed in liquid lithium at two different 
temperatures: T — 475 K (slightly above the melting point) and 600 K. Taking advantage of the 
absence of any kinematical restriction and incoherent contribution, and pushing the instrumental 
resolution up to 1.5 meV, it was possible to perform an accurate investigation of the dynamic 
structure factor S{Q, E) in the wavevector range from 1 to 110 nm~^. For Q smaller than Qm — 25 
nm~^ , the position of the main peak of the static structure factor, a detailed analysis of the lineshapes 
shows that any picture of the relaxation mechanisms based on a simple viscoelastic model must be 
abandoned. All the spectral features can instead satisfactorily be accounted for by including both 
fast and slow relaxation processes. The physical origin of the slow relaxation is associated to the 
structural rearrangement, while the local nature of the fast one is extensively discussed. At larger 
Q values a gradual crossover from the strongly correlated to single particle dynamics occurs, with 
an important weight provided by quantum effects. 



PACS numbers; 61.25.Mv, 61.20.Lc, 67.40.Fd 
I. INTRODUCTION 

The dynamics of liquid metals has been extensively in- 
vestigated in the recent past with the main purpose of 
ascertaining the role of the mechanisms underlying both 
collective and single-particle motions at the microscopic 
level. In the special case of collective density fluctua- 
tions, it is known that well-defined oscillatory modes can 
be supported even outside the strict hydrodynamic re- 
gion. In molten alkali metals, this feature is found to 
persist down to wavelengths of one or two interparticle 
distances, making these systems excellent candidates to 
test the various theoretical approaches developed so far 
for the microdynamics of the liquid state. As a conse- 
quence, since the pioneering inelastic neutron scattering 
(INS) study by Copley and Rowe jl] in liquid rubidium 
the interest in performing more and more accurate exper- 
iments is continuously renewed: INS investigations have 
been devoted to liquid cesium , sodium |^ , lithium , 
potassium 1^ and again rubidium 

Up to a few years ago the only experimental probes 
adequate to access the interparticle distance region in 
collective dynamics were thermal neutrons. With this 
technique foundamental results have been achieved in the 
field of condensed matter structure and dynamics. Un- 
fortunately, in several systems their use for the determi- 
nation of the dynamic structure factor S{Q^ E) becomes 
extremely difficult (if not impossible) for two reasons. 
The first one reflects the presence of an incoherent con- 
tribution to the total neutron scattering cross section. In 
liquid sodium, for example, the incoherent cross section 
dominates; even in more favorable cases (Li, K) at small 
Q the intensity of the collective contribution is low, and 



its extraction requires a detailed knowledge of the sin- 
gle particle dynamics. The second reason is dictated by 
the need of satisfying both the energy and momentum 
conservation laws which define the (Q — E) region ac- 
cessible to the probe 0. Roughly speaking, when the 
sound speed of the system exceeds the velocity of the 
probing neutrons (~ 1500 m/s for thermal neutrons) col- 
lective excitations cannot be detected for Q values below 
Qmi the position of the first sharp diffraction peak of the 
sample, namely just in the most significant region for the 
collective properties. 

By virtue of the m~^/^ dependence of the isothermal 



sound speed cq (see section HA) the higher the atomic 
number of the system, the wider is the kinematic re- 
gion accessible to neutrons, so that accurate INS data 
are available for rubidium (cq ~ 1260 m/s) and ce- 
sium (co - 970 m/s) §, while more difficulties are met 
in the case of lighter metals. In particular, lithium rep- 
resents the most critical case due to its high sound speed 
(co ^ 4500 m/s) and to the equality of the coherent and 
incoherent cross sections: for this reason INS aiming to 
the study of collective properties of Li represents a very 
hard challenge jsj . From a general point of view the main 
outcome of all these INS experiments, as far as collective 
properties are concerned, is the evidence of inelastic exci- 
tations in S{Q, E) that, in the specific cases of rubidium 
and cesium, exhibit a positive dispersion, i.e. an increase 
of LOmiQ) (the position of the S{Q,huj) peaks) with re- 
spect to the values implied by the hydrodynamic sound 
speed. About the details of the relaxations driving such 
dispersion not much can be inferred at this level: the ex- 
perimental results have been necessarily analyzed within 
simple models such as the damped harmonic oscillator 
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1^ or the kinetic model H,^, suitable to extract reliable 
and resolution-corrected information on the peak posi- 
tions but not about the lineshape features. Some extra 
efforts have been done in the case of cesium where in- 
formations about an average relaxation time have been 
extracted 

A very useful tool, complementary to the "tradi- 
tional" experimental facilities is the numerical simulation 
technique, in particular Molecular Dynamics (MD): the 
choice of a realistic interatomic potential - i.e. a potential 
model able to reproduce structural properties - allows the 
determination of the dynamics of the system via the inte- 
gration of the classical Newton equations. In this frame- 
work the single particle and the collective properties can 
easily be investigated within technical restrictions due 
to the finite box size (defining the minimum accessible 
wavevector) and computation time (related to the statis- 
tical quality and to the energy resolution of the calculated 
spectra). Broadly speaking the features of the atomic col- 
lective motion i.e. the details of the S{Q, E) lineshape, as 
outcome of MD run, turns out to be less noisy and more 
straightforward than the correspondent INS results: no 
absolute normalizations are required, no mixing between 
coherent/incoherent dynamics occurs and, above all, ba- 
sically no resolution corrections are needed. 

After the first experiments on rubidium, a considerable 
number of MD works have been published in this field, 
giving a valuable support to the experimental measure- 
ments . On the specific case of liquid lithium, the 
performed MD simulations show features common 
to all the other molten alkali metals, like the presence of 
positive dispersion - a particularly relevant result that 
could have not be achieved with neutrons for the previ- 
ously mentioned reasons. At temperature around T™ the 
increase of the sound velocity is of the order of 20%, and 
not much is known about its microscopic origin. 

In the recent past, the development of new synchrotron 
radiation facilities opened the possibility of using X-rays 
to measure the S{Q, oj) (which is proportional to the scat- 
tered intensity) in the non-hydrodynamic region; in this 
case the photon speed is obviously much larger than the 
excitations velocity, and no kinematic restriction occurs. 
Moreover, in a monatomic system as lithium, the scat- 
tering cross section is purely coherent and so it is directly 
associated with the dynamic structure factor. Some ex- 
periments have been performed on liquid lithium |l4|-[l^] 
with progressively increasing resolution; it has been pos- 
sible to show the existence of propagating collective ex- 
citations, but due to resolution limitations (never below 
AEpwHM ^ 12 meV), no detailed information about 
the hneshape could be extracted. 

In this work we report the results of Inelastic X- 
ray Scattering (IXS) experiments performed on liquid 
lithium with very high energy resolution, adequate to 
probe accurately the detailed features of S{Q,E) - this 
allows to investigate whether the lineshape dependence 



on momentum transfer can be interpreted within the phe- 
nomenology of one or more relaxation processes. In par- 
ticular in Section II we review the basic theoretical frame- 
work adopted for the interpretation of the IXS results. 
The latter are reported in Section III together with a 
brief account of the experimental setup. Section IV is 
devoted to the data analysis and discussion. The main 
outcome of the present paper is finally summarized in 
Section V along with some concluding remarks. 

II. BASIC THEORY 

A. Definitions 

In the classical limit, the basic time correlation probing 
collective dynamics in a monatomic fluid (iV particles 
with mass m) is the intermediate scattering function 

F{Q,t) = {l/N)J2{e''^"'''^"^^'^"''^'^) (1) 

where ri{t) denotes the position of the j-th particle at 
time t. The dynamic structure factor S{Q,lu) is the fre- 
quency spectrum of F{Q, t), while the structural features 
are accounted for the initial value F{Q,t = 0) = S{Q). 
The quantity F{Q, t) then obeys the equation 

F iQ,t)+LoliQ)F{Q,t) (2) 

+ [ M{Q,t-t') F {Q,t')dt' = 

where 

u;2(Q) = KBTQ^mSiQ) (3) 

is the second classical normalized frequency moment of 
S{Q,uj), while M{Q,t) is the so called memory function 
of the system, related in some way to the details of the 
Hamihonian In the hmit Q ^ 0, uj^{Q) c^Q^, 

where cq is the isothermal sound velocity; hence the 
quantity cq (Q) = iOo{Q)/Q can be interpreted as a suit- 
able generalization of cq to finite wavevectors. 

From Eq.(H) a formally exact representation of the 
Laplace transform of F{Q, t) can be written as ||7|j20|| 

/>OC 

F(Q,z)= / dte-''F{Q,t) (4) 
Jo 

\ lz+M(Q,z)]j 

From the knowledge of F{Q, z) one straightforwardly 
obtains S{Q,uj) — (l/7r)i?eF(Q, z — iuj) in terms of 
the real (Af) and imaginary {M") parts of the Fourier- 
Laplace transform of the memory function: 
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S{Q,u;) 



S{Q) 



Lu^o{Q)M'{Q,Lo) 



(5) 



The spectral features of the dynamic structure factor 
can be characterized by its frequency moments 17^"'' (Q) = 
/ uj^S{Q,uj)duj, where, for a classical system, only the 

even frequency moments (such as rt^g\Q) = S{Q) and 
n'-^^Q) = S{Q)uj^{Q)) are different from zero. 

In the following, we shall also find convenient to con- 
sider the "longitudinal current spectrum" defined as 
Cl(Q,(^) = (cj^/Q^) S{Q,uj). The presence of the factor 
uj'^ wipes out the low frequency portion of the dynamic 
structure factor, and consequently emphasizes the gen- 
uine inelastic features of S'(Q,w). After its definition 
and Eq.(Q), it is readily seen that the Laplace transform 
Cl{Q, z) satisfies 



Cl{Q,z) ^ -z[zF{Q,z) ~ S{Q)] 

= [z+[loI{Q)/ z] + M(Q,z)}"' 



(6) 



Again, the spectrum Cl{Q,uj) can be expressed as 
{l/Tr)ReC'L{Q, z = ilu). Then the position and the width 
of the inelastic peaks in Cl((3, lu) are determined by the 
poles of Cl{Q, z). 

The above definitions are valid for a classical system. 
In the case of our interest, the main effect of quantum- 
mechanical corrections stems from the well-known in- 
equality of the positive and negative-frequency parts of 
the spectra, connected by the detailed balance factor 
e^^^ . Additional sources of non-classical behavior, such 
as those associated with a finite value of the dcBocr wavc- 

1/2 

length A = {2nh^ / mK bT^ , are small (in the explored 
lithium states A is only 0.11 times the average interpar- 
ticle distance) and can safely be neglected. Since the 
effects of the detailed balance are clearly visible in the 
experimental IXS spectra, we briefly discuss a possible 
procedure to account for this quantum feature in a con- 
sistent way, while preserving the inherent advantages of 
the classical description. In doing this, for the sake of 
clarity we shall denote all the previous classical quanti- 
ties with the subscript c/, while the notation q will refer 
to the quantum case. 

The natural theoretical counterpart of the classical 
density correlation function is the so called Kubo canon- 
ical relaxation function 



Kg{Q,t) ^ -^y2 I dX ( e-'^-'-^'^^'e--^'' e'^-'-^'-'-'e^ 



/3 



,-jQ-?.(0)„-AH iQ-¥j(t) AH 



(7) 



where (3 ~ 1/KbT and the angular brackets denote 
a quantum statistical average. In the classical limit 



{(3^0, h ^ 0) the operators A become classical com- 
muting dynamical variables and Kq(Q,t) —f Fci{Q,t). 
It can be shown that Kq{Q,t) is a real even func- 
tion of time, so that its spectrum Kg(Q,uj) is an even 
function of frequency. On the other hand, the ex- 
perimental scattering cross section involves the Fourier 
transform Sq{Q,Lu) of the quantum density correlator 
Fq{Q,t) = {l/N)J2i,j {e~'^''^'^°'>e^^-'^^^^. The relation 
between Sq{Q,u)) and Kq{Q,uj) reads 

and satisfies the detailed balance condition. Moreover 
the condition 



-1) 



(8) 



relates the even frequency moments of Kq with the odd 
ones of Sq. The same memory function framework of Eq. 
(^ can be phrased for the Kubo relaxation function and 
for its Laplace transform Kq{Q, z). 

By virtue of all these properties, in a situation where 
the quantum aspects not associated with detailed balance 
are marginal, it is reasonable (although not strictly rigor- 
ous) to identify the spectrum Kq{Q,uj) with the classical 
quantity Sci{Q,uj) so that 



(3huj 
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Sci{Q,uj) 



(9) 



Having established such a correspondence, from now 
on we will drop out the subscript cl and refer to the 
classical quantities as in fact done at the beginning of 
this section. 



B. The memory function features 

The memory function M{Q, t) accounts for all the re- 
laxation mechanisms affecting collective dynamics, and 
consequently is the central quantity in most theoretical 
approaches. From straightforward algebra, the initial 
value of M{Q,t) is related to the spectral moments of 
S{Q,uj) by: 



M(g,i = 0) 



4^)(Q) 



(10) 



An exact expression of fl^g\Q) exists, but involves both 
the derivatives of the interparticle potential and the pair 
distribution function j?). It is usual to define ijj\{Q) = 
nf{Q)/^lf{Q) and A2(g) = loHQ) ~ uj^Q), their 
meaning being evidentjrom the following argument. For 
sufficiently large \z\, M{Q,z) ~ M{Q ,t = 0)/z and E q. 
(P) is seen to have poles at z = ±iy\v^{Q) + A^{Q) — 
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±iu}L{Q), showing that the frequency lul{Q) character- 
izes the instantaneous collective response of the liquid 
at the wavevector Q. Similar remarks can be made 
for the generalized infinite- frequency velocity Coo{Q) = 

>^l{Q)/Q- 

In liquid systems, in the long time limit, one expects 
that M{Q,t oo) approach zero value, therefore, re- 
gardless the details of its shape, M{Q,t) it is expected to 
decay over a certain timescale t{Q) that, for the sake of 
simplicity, can be defined as t{Q) = M(Q)/A^((5), being 
= M{Q,t)dt. It can be of considerable inter- 
est to point out the asymptotic behaviors of Eq. (^) in 
the opposite regimes t{Q)llJo{Q) << or >> 1. In the first 
limit (M (Q, t) « 2M{Q)S{t)) Eq. (|) reduces to the equa- 
tion for a damped harmonic oscillator, and it can be easily 
proved that M'{Q,uj) = M{Q) and M"(Q,w) = 0, which 
means two inelastic peaks in the (in the current 

spectra they are centered at ±wo(Q)) damped with a 
factor M{Q). In the opposite limit, t{Q)u!o{Q) >> 1, the 
decay of F{Q,t) is much faster than t{Q) and M{Q,t) 
appearing in the convolution integral of Eq. (||) can be 
considered constant, M{Q,t) « A'^{Q). Therefore Eq. 
(B) becomes non homogeneous and the solution of Eq. 
(5) reduces to a harmonic oscillator of frequency ujl{Q) 
with no damping, plus a sharp elastic line. In this ex- 
treme case, the intensity ratio between the elastic and 
inelastic lines (the Debye Waller factor) in S{Q,uj) is 
/(Q)-l-c.2(Q)K(Q). 

The true time dependence of the memory function is a 
priori unknown, and from the very start it is convenient 
to separate in M{Q,t) the decay channels which explic- 
itly involve couplings to thermal fluctuations ( Mth{Q, t) 
) from those directly associated with longitudinal density 
modes ( ML{Q,t) ). A convenient way to perform this 
splittingis by the use of the generalized hydrodynamics 
theory [0,0. In the simplest version of this approach, 
one writes 

M{Q,t) = ML{Q,t)+MthiQ,t) (11) 
= Al{Q)mL{Q,t) + AUQ)mth{Q,t) 

with Ai(g) ^ uliQ) - liQHiQ), AliQ) = [7(Q) - 
1]ujq{Q). Here 'y{Q) is a generalization of the specific 
heat ratio 7 — Cp/Cy to finite wavevectors. The relax- 
ation processes underlying the dynamics are accounted 
for by the normalized quantities mL{Q, t) and mth{Q, t), 
defined in such a way that mL{Q, 0) = mth{Q, 0) = 1. A 
straightforward generalization of ordinary hydrodynam- 
ics suggests for the thermal contribution the following 
form 

mth{Q, t) « exp[-a{Q)Qh] (12) 

where a{Q) can be viewed as a finite Q generalization of 
the quantity Dt = n/nCy, being k, the thermal conduc- 
tivity. The role of the thermal contribution in the dy- 
namic structure factor of liquid metals will be discussed 
in the following section. 



In contrast with the thermal decay channel, no guid- 
ance for Mi((5,t) is provided by ordinary hydrodynam- 
ics. In the latter, one implicitly assumes that 



Ml{Q ^ 0,t) » 2{T^Llnm)Q''5{t) 



(13) 



where the longitudinal viscosity coefficient 77L is related 
to the ordinary shear and bulk viscosities by t/l = 
(4/3)?7 -I- rjB- Clearly in the limit Q ^ the ratio 
Ml{Q, z = 0)/Q^ approaches rjL/nm. 

The simplest way to go beyond the hydrodynamic re- 
sult (|l^ ) is to allow for a finite decay rate of ML{Q,t): 



ML{Q,t)^Al{Q)e 



-t/r(Q) 



(14) 



where it has been additionally assumed an exponential 
lineshape for the memory function. Although this has 
the advantage of analytical simplicity when dealing with 
Fourier transform, a drawback of this ansatz lies in the 
violation of some basic short-time features of the memory 
function (such as a zero derivative at t = 0) , causing the 
divergency of Jl^"'' for n > 6. 

Eq. (|l4|) yields the so-called viscoelastic model for 



Since as Q ^ mL{Q,z = 0)/Q can 
— Cq]t{Q 0), the requirement that this 



SiQ,u;) 
written as 

coincides with rj^/nm shows that the time t{Q) must be 
finite as Q ^ 0. Such a connection with viscous effects 
justifies the physical interpretation of the rate 1/t{Q) 
as a parameter giving an overall account of all relaxation 
processes by which the longitudinal response of the liquid 
is affected by time-dependent disturbances. In particu- 
lar, for slow perturbations developing over a timescale 
t » ^iQ) the system can adjust itself to attain local 
equilibrium and the usual viscous behavior. In contrast, 
for disturbances fast enough that t ^ t{Q) the hquid 
responds instantaneously, with a solid-like (elastic) be- 
havior. The crossover between these limiting situations 
(times t « t{Q), or frequencies uj such that ujt{Q) « 1) 
is ultimately responsible for the gradual changes often 
detected in the sound dispersion of several liquids at in- 
creasing wavevectors. 

Similar considerations can be applied to the Mth{Q, t) 
contribution which involves the timescale Tth — 
Dt^{Q)Q^; in this case the "elastic response " is achieved 
in the hydrodynamic regime: lOoiQ 0)r((5 ^ 0) a 
co/Q 00, while the "viscous response " emerges at 
short wavelengths. The crossover between the two lim- 
its embodies the transition between the adiabatic and 
isothermal regimes, and refiects the actual possibility of 
density fluctuations to decay by the thermal channel: at 
high frequencies there is no time for such conversion, 
and the fluctuations evolve without heat transfer (higher 
sound speed, no attenuation) while in the low frequency 
limit there is enough time for the system to equilibrate 
so that an isothermal dynamics takes place (lower sound 
speed, damped excitation). 
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These quaUtative considerations are supported by the 
evolution of the poles of Cl(Q, z) at increasing wavevec- 
tors. Eqs. (|ll|), (ph and in fact imply that 



Mi(Q,i)=A2(g)e-*/%+A2(Q)e 



(17) 



M{Q,z) 



z + 1/t{Q) ' z + a{Q)Q^ 



A?,(g) 



(15) 



In the hydrodynamic regime {Q — > 0) both the in- 
equalities l^l ^ l/'''(0) and \z\ » a{Q)Q^ are satis- 
fied, and the poles are approximately located at z « 
±i^/Lu^{Q ^0) + A2^(g 0) = ±i Cs Q, where Cs = 
Co is the adiabatic sound velocity. At larger wavevec- 
tors different situations may arise: the one of our interest 
is \z\ <^ a{Q)Q'^ , which is readily seen to cause a shift 
of the imaginary part of the poles toward an isothermal 
behavior, i.e. uj ~ ±wo((5). The last limiting case occurs 
whenever the frequency w becomes distinctly larger than 
[1/t{Q)], which eventually yields poles at 2; = zLujl{Q), 
as already noticed. Although appealing, the simplicity 
of the viscoelastic model can be deceptive. First of all, 
the model itself provides no clue for the physical origin of 
the decay mechanisms leading to the rate 1/t{Q). Sec- 
ondly, even ignoring this aspect and treating 1/t{Q) as 
a fitting parameter, the practical results of the model 
are rather unsatisfactory (see Fig. |5| in the following). 
These drawbacks of the simple viscoelastic model have 
been theoretically predicted, and somehow reported in 
a number of liquids through a detailed analysis of MD 
spectra [0j23|. However, to our knowledge, no cutting 
edge analysis exists in the case of experimental measure- 
ments of the coherent dynamic structure factor, proba- 
bly because of the afore-mentioned limits of the neutron 
technique. 

The obvious remedy is to modify the simple ansatz 
(0) by allowing a more sophisticated decay of Ml{Q, t). 
We adopted the following two-exponential ansatz: 



Mi(g,i) = Ai(Q) 



(l-a(g))e-^i^'^^*+a(g)e-^^^'5)*' 

(16) 



where the rate 71(g) is chosen to be larger than 72(g), so 
that the dimensionless factor a{Q) measures the relative 
weight of the "slow" decay channel. Besides being more 
fiexible than the viscoelastic model, we shall see that the 
ansatz (^|) has the much more important merit that the 
presence of two different timescale does in fact have a 
definite physical interpretation. 

The previous ansatz has been proposed phenomenolog- 
ically several years ago and tested against the MD results 
in LJ fluids |^^, it is interesting to note that an expres- 
sion analogous to Eq. ( p!6| ) has been implicitly introduced 
in the viscoelastic analysis of Brillouin Light Scattering 
(ELS) spectra of glass forming materials (see for example 
ref. j2^]). In fact, in these works, the general expression 
of ML{Q,t) for a two times exponential decay is always 
expressed as 



with explicit reference to the so-called a- (structural) re- 
laxation process as responsible for the long lasting tail, 
and to the /i— microscopic process as additional, faster, 
relaxation dominant over a very short timescale. To be 
more precise, in the ELS spectral window the condition 
LUT^ << 1 holds, so that the approximation 

Mi(g,t) = 2A2(Q)r^J(t)+A2 (g)e-*/-» 

is customarily adopted as a fitting model. On the con- 
trary, we will show in this work how such an approxi- 
mation is no longer tenable in the case of liquid Lithium 
at the IXS frequencies. The origin, at the atomic level, 
of this fast decay channel is still an open issue: one of 
the purposes of the present paper is to point out the fea- 
tures of such relaxation and to inspect its nature. The 
rapidly decaying portion of Afi(g,t) is customarily at- 
tributed to largely uncorrelated collisional events, sim- 
ilar to those occurring in a dilute fluid. In addition, 
at the high densities typical of the liquid state, non- 
negligible correlations among the collisions can be ex- 
pected, making no longer valid an interpretation only 
in terms of "binary" collisions. Albeit the magnitude 
of the correlation effects is relatively small and their 
buildup slow, once established their decay is even slower, 
and for t > l/7i(g) = this relaxation channel dom- 
inates the decay of ML{Q,t), which consequently may 
exhibit a small but long-lasting " tail" . The ansatz ( |T^ ) 
can incorporate most of this physics: on the basis of 
the latter, one may reasonably anticipate that a{Q) ^1, 
and that the time l/72(g) = is distinctly longer than 
1/71 (g) = fi- I'^ t^i^ picture, the best fltted values of 
the viscoelastic rate 1/t{Q) clearly represent some sort 
of "weighted average" between 71(g) and 72(g). Fi- 
nally, we may argue that at increasing Q (namely, over 
a shrinking length scale) the magnitude a{Q) of correla- 
tion effects should decrease, and that at higher temper- 
atures the value of a{Q) at a given wavevector should 
equally decrease. On a general basis, the requirement 
that liniQ^Q M]^{Q, z = 0)/Q^ — > rjL/nm now takes the 
form 



(c 



cl) 



(l-Q(g^O)) ajQ ^ 0) 
7i(g^0) 72(g^0) 



rj^/nm 
(18) 



The refined model (16) yields a dynamic structure factor 
given by 



xRe 



S{Q) 



(19) 



^o'(g) 



AL(Q) 



ii^+'ri{Q) i'^+l2{Q) iuj+a{Q)Q'^ 
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III. THE EXPERIMENTS 

In this work we report the determination by IXS of 
the dynamic structure factor of hquid hthium in the 
wavevector range from 1.4 to 110 nm^^ , corresponding 
to Q/Qm ~ 5 • 10~^ 5. The scanned energy range has 
been settled at each Q in order to detect all the scattered 
signal up to the tails region, where it becomes compara- 
ble to the background. 

The experiment has been performed at the very high 
energy resolution beamline ID16-BL21 at the European 
Synchrotron Radiation Facility. A monochromatic beam 
of 10^ photons/s is obtained from a cryogenically cooled 
Si(lll) double crystal followed by a high energy resolu- 
tion monochromator operating in back-scattering geom- 
etry at selectable Bragg reflections The scattered 
beam is collected by perfect spherical silicon crystal an- 
alyzers operating in back-scattering and Rowland circle 
geometry at the same reflection order of the monochro- 
mator. They were obtained by gluing ~ 12000 perfect 
crystals of 0.6 x 0.6 x 2 mm'^ on a spherical blank [^ . 
The overall energy resolution has been measured using 
a plexiglass scatterer at a Q value corresponding to the 
maximum of its structure factor where the diffusion is 
dominated by the elastic contribution, and it turns out 
to have approximately a Lorenzian shape. 

Different configurations have been adopted depending 
on the explored Q-range. Below Qm (25 nm~^) we used 
simultaneously five-analyzers mounted on a horizontally 
rotating arm and with a fixed offset of 1.5°. Either the 
Si(9 9 9) or Si(ll 11 11) reflections, corresponding to 
energy resolutions of 3.0 meV and 1.5 meV (FWHM) re- 
spectively, have been selected. At higher Q values, where 
the spectra start to broaden, we used a vertical scattering 
geometry based on a single analyzer housed on a rotat- 
ing arm that can reach more than 150° scattering angle. 
In this configuration we selected the Si (7 7 7) refiection 
with an overall energy resolution of 8.5 meV FWHM. The 
Q resolution was ~ ±0.2 nm~^ and ~ ±0.35 nm^^ in the 
horizontal and vertical geometries respectively. The mo- 
mentum transfer is related to the scattering angle by the 
relation Q = 2kosin{9s/'2). 

Energy scans have been done by varying the temper- 
ature of the monochromator with respect to that of the 
analyzer crystals. The absolute energy calibration be- 
tween successive scans is better than 1 meV. Each scan 
took about 3 hours, and each Q-point spectrum has been 
obtained from the average of 2 to 8 scans depending on 
the Bragg-reflection order and on the Q-transfer. The 
data have been normalized to the intensity of the inci- 
dent beam. In the Q—E region of interest, empty vacuum 
chamber measurements gave either the flat electronic de- 
tector background of 0.6 counts/min or, in the angular 
region corresponding to 9 < Q < 13 nm^^, a small elas- 
tic line due to spurious reflection of the beam entering 



the analyzers. These background signals have been sub- 
tracted from the data before any other data treatement. 
The chamber kapton windows (each 50 /xm thick) gave 
no detectable contribution to the scattered flux. 

The liquid lithium uncapped container was made out 
of austenitic stainless steel in contact with a resistive 
heater, used to kept the liquid at constant temperature. 
We selected two different values: one slightly above the 
melting point, T = 475 K, and the second at T = 600 K; 
we have not been able to reach higher temperatures with 
our experimental setup due to the lack of windows. The 
20 mm long sample, kept together by surface tension, was 
maintained in a 10~^ bar vacuum. The lithium has been 
loaded in an argon glove box. 

As far as experimental aspects are concerned, the en- 
ergy transfer E between the probe (X-rays) and the sam- 
ple (excitations) has always been measured in meV. To 
make contact with the traditional notations adopted in 
computer simulations we will refer from now on only to 
the frequency uj = 27ru = E/h expressed in ps~^ (a nu- 
merical factor of « 0.659 ps~^ /mcV accounts for such a 
conversion) . 

The IXS spectra of liquid lithium below Qm (horizon- 
tal geometry) are reported in Fig. |^ at the indicated 
Q values. The low Q data show the Brillouin triplet 
structure with the energy of the inelastic peaks increas- 
ing with Q up to a Q value of 12 nm^^. This value is close 
to Qml'2, as deduced from the Siff) reported in Fig. ||. 
The dispersion up to Qm/'^ can therefore be interpreted 
as that of the longitudinal acoustic branch in a pseudo 
first Brillouin zone (BZ). Furthermore, similarly to what 
is found in the second BZ of a crystal, we observe that, 
also in liquid lithium, the energy of the acoustic modes 
decreases with increasing Q from Qm/'^ to Qm- Increas- 
ing Q above Q™, i. e. in the "third" or higher BZs, the 
spectrum increasingly broaden and distinct peaks are no 
longer observable. At the highest Q- values one finds that 
the dynamic structure factor becomes a symmetric peak 
centered at frequencies larger than w = 0, the behav- 
ior expected for a quantum free particle (in the classical 
limit no recoil energy is expected). Such evolution can 
be clearly followed in Fig. Beside the observation of 
dispersion in a first and second pseudo- BZs, it is also im- 
portant to note that the broadening of the excitations 
increases monotonically with Q, to the extent that near 
the end of the second pseudo-BZ a well defined inelastic 
peak is no longer visible. 

In order to extract quantitative information from the 
experimental intensity, i.e. to perform measurements of 
Sq{Q, Lu) on an absolute scale, the most direct way would 
be to use a reference scatterer as customarily done in neu- 
tron experiments. In IXS such a procedure is extremely 
difficult because of the Q-dependence of the form factor 
and, in our case, for the different efficiencies of the ana- 
lyzers. For these reasons we preferred to use an alterna- 
tive indirect method based on the knowledge of the sum 
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rules of Sq{Q, w): in particular for the first two frequency 
moments we have 

4"' = / S,iQ,Lu)diu = S{Q), 

n'-g^ = J LuSq{Q,uj)duj = hQ^/2m. 

where the second equality follows from Eq. applied 
for n = I. The measured raw intensity is related to the 
dynamic structure factor through 

I{Q,Lu) = A{Q) J duj'Sg{Q,Lu')R{u;~Lj') (20) 

where R{lu) is the experimental resolution function and 
A{Q) is a factor taking into account the scattering ge- 
ometries, the experimental setup and the lithium atomic 
form factor. The first moments of the experimental data, 
flf'^ and fl'^P , and those of the resolution function, Sl^-* 
and Sl^-* , are related to fl^g^ and fl^g'^ by: 

nf^ = A{Q)nfnf, 

From the previous equation one derives that 

= (21) 

This procedure has been adopted to establish an ab- 
solute scale for Sq{Q,Lu) using the experimentally deter- 
mined I{Q,u}) and R{uj). Its reliability is shown in Fig. 
^ where we obtain an excellent agreement between the 
Sq{Q) values obtained by Eq. ( pl| ) and those derived by 
a MD simulation at T = 475 K |2|,||l. The reported 
MD data are the outcame of a simulation made using 
the Price- Tosi interaction pseudopotential already tested 
against neutron diffraction data (see ref for further 
details). 

IV. DATA ANALYSIS 

Lithium has been chosen because, among the simple 
monatomic liquids, is the one that is best suited to IXS. 
Indeed, its low mass gives recoil energies observable in the 
considered Q-range, and its low atomic number and large 
sound velocity give optimal signal within the available en- 
ergy resolution, compensating for the large form factor 
decrease at high Q values. The I{Q,uj) spectra, reported 
to their absolute scale exploiting the zeroth and first mo- 
ments sum rules, show the transition from a triplet to 
a Gaussian. The maxima, ujm{Q), of the longitudinal 
current spectra {uj'^S{Q,uj)/Q'^) show an almost linear 
dispersion relation at low Q - typical of a sound wave 
- and a completely different dependence in the high Q 



limit, where they approach the parabolic dispersion of 
the quantum free particle. Between these two regions, 
^m{Q) exhibits oscillations which are in phase with the 
structural correlations as observed in S{Q). 

The low Q region is of particular interest for its large 
number of informations. The so called positive disper- 
sion, a fingerprint of the relaxation dynamics of disor- 
dered systems, either hquid |l|,^,|§gl or glassy |^,|ll: 
reflects the host of mechanisms driving such dynamics, 
their timescale and, at a careful inspection, their nature. 
We will start our analysis from here. 

A. The "low Q region" 

Looking at the spectra in Figs. |l],|[ as also confirmed 
by the results reported in this section, the memory func- 
tion approach can profitably be used approximately up 
to Qm. We adopted the following procedure. The classi- 
cal result of Eq. (||) has been modified according to Eq. 
(^) to take into account of the detailed balance effects. 
The resulting Sq{Q,Lu) has then been convoluted with 
the instrumental resolution according to ([20|). The 
function arising from the difference between the outcom- 
ing function and the experimental data has finally been 
minimized by a standard Levenberg-Marquardt routine. 
Among the fitting parameters, S{Q) and i^o{Q) are inde- 
pendent of the specific model for the memory function, 
and their value is basically known: the first has been cal- 
culated from the first moment sum rule as discussed in 
Sec. Ill (we have already discussed the approximation 
Sq{Q) = S{Q)), so that it has been kept fixed during 
the iterations. Then the second quantity ujq{Q) is simply 
deduced from its definition (see Eq. (^)). It is worth 
to point out that uJq{Q) cannot be taken as the second 
frequency moment of I{Q,lu) for two main reasons. The 
first one, related to the effects of the resolution, may be 
in principle overcome by a procedure similar to the one 
leading to the determination of S{Q), but with the fur- 
ther difficulties arising from the almost-Lorcnzian shape 

(2) 

of the resolution that leads to a diverging ri)j . Most 
severe is the non-invariance of the second moment under 
the transformation (|^). For the sake of clarity we report 
in Fig. ^ the theoretical value u!q{Q) compared to the 
second normalized moment of the experimental I{Q,uj) 
at T = 475 K: the latter is always larger (this may be 
due also to the small but non-vanishing background in 
the experimental data) and the agreement gets worse at 
decreasing wavevectors, when the resolution effects be- 
come dominant. 

All the other fitting parameters obviously depend on 
the specific model for M{Q,t). Let us firstly exam- 
ine the spectral relevance in liquid Li of the coupling 
to thermal fluctuations, accounted for by the contribu- 
tion MfhiQ, t) to the total memory function. At not too 
high temperatures in the molten phase, all liquid alkali 
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metals are characterized by a specific heat ratio 7 « 1, 
the typical values giving 7 — 1 ~ 0.1. Outside the hy- 
drodynamic (Q — s- 0) regime (and, in particular, in the 
wavevector range probed in this work) there are several 
indications that the difference 7(Q) — 1 is even smaller 
[^,|2|. Despite its modest strength in the present Q- 
range, the role of thermal relaxation in liquid metals 
is strikingly different from the one played in ordinary 
non-conducting liquids. A simple estimate of the de- 
cay rate embodied in the Mth{Q,t) reported in Eq. ( p^ ) 
based on the experimental values of the thermal conduc- 
tivity and of the specific heat shows that, in liq- 
uid lithium, uj <C {K/nCv)Q^ for most of the spectral 
range of interest here, and in particular for the range 
characterized by the presence of inelastic peaks (see Fig. 
|l|). This indicates that in liquid Li (but a similar ar- 
gument is expected to hold in general for molten met- 
als) the afore-mentioned transition from the adiabatic to 
the isothermal response occurs at wavevectors well he- 
low those probed in the present experiments. Indeed, for 
wavevector-independent thermal parameters (a reason- 
able assumption at so low Q) we find for the threshold 
wavevector uJo{Q*)Tth{Q*) = 1 a value of Q* ~ 0.2 nm^^ 
at T = 475 K, varying only of few percents at T = 600 K. 
Of course, this does not mean that the presence of ther- 
mal fluctuations is not influent: as shown for example by 
Eq. ( pi), an d as can be inferred by the simple arguments 
of Sec. [I A, they may still give some contribution to the 



real part of the memory function, i.e. ultimately to Bril- 
louin peaks widths. Such linewidth contribution gives 
Sujth = ^Ih^h = (7(Q) - 1)co{Q)/Dt{Q) ^ 0.056 ps-i 
(0.07 ps"i) for T = 475 K (T = 600 K) (in the Q 
limit) and is not explicitly dependent on Q. As already 
remarked, such contribution may be of some relevance 
only at the first two or three Q values investigated here, 
but the strongly Q-dependent broadening of the inelastic 
components has to be ascribed to some other dynamical 
process. On the basis of the same argument, it is worth 
pointing out that the hydrodynamic value of the sound 
speed is not yet reached at the lower Q values of this 
work: as far as the positive dispersion is concerned, in 
the present case, the lower edge of the transition is the 
isothermal rather than the adiabatic value. As previously 
noted, a completely different scenario characterizes in- 
stead non- metallic fluids: due to the much lower thermal 
diffusivity, the opposite thermal regime is experienced all 
over the collective excitation region and the final low Q 
result is an adiabatic response without any damping of 
the sound waves and with a narrow elastic line. 

Having ascertained the role of thermal effects, we now 
come to the more important implications of the fitting 
analysis for the dynamics of mL{Q,t). In the following 
we will refer to Fig. || where the specific situation of the 
spectrum at Q = 7 nm~^ has been chosen to enlighten 
the lineshape details associated with the different contri- 
butions to the memory function. 



As a zero-level approximation we consider the line- 
shape obtained by the simple hydrodynamic expression 
which only accounts for thermal decaying processes, em- 
bodying all the other relaxation channels as instanta- 
neous: beside the use of Eqs. (|ll|,|l^, this is simply 
done assuming that mL{Q,t) cc S{t). As in our case 
even the thermal relaxation is quite rapid on the X-rays 
timescale, the net result is a lineshape close to the one of 
the damped harmonic oscillator model. It is therefore not 
surprising (see Fig. ^ that the agreement with experi- 
mental data is very poor (no elastic peak is reproduced) . 
Consequently, the assumption of an instantaneous decay 
of mL{Q,t) is unten able. 

Following section [IB , as a next step, we generalize 
the memory function by the expression (|l5| ) i.e. the vis- 
coelastic model. All the merits and the drawbacks of 
this model are again illustrated in Fig. ^. Leaving the 
time t{Q) and the strength A|(Q) as fitting parameters, 
it is possible to obtain a rather good agreement for the 
position of the inelastic peaks, and consequently for the 
main features of the dispersion of sound-like excitations. 
However, one notes clear discrepancies in the quasi- elas- 
tic region at small lu. The fitting procedure yields with 
comparable values two different sets of parameters, 
corresponding to two different timescale driving the dy- 
namics of the system. More precisely, in one case the 
broader part of the quasielastic and the Brillouin line- 
shape are well reproduced but the narrow quasi-elastic is 
missing, while in the other one the narrow elastic peak 
is well fitted but the Brillouin lines are too sharp. In 
other words, the viscoelastic model cannot account for the 
double-sloped shape of the quasielastic region. Although 
in Fig. ^ we refer to Q = 7 nm^^, the situation at most 
wavevectors is similar. 

The inescapable conclusion is that there are at least 
two different timescale, each giving a contribution to the 
width of the "quasi-elastic" region, and accounting at the 
same time for the inelastic lineshape. On these grounds, 
the more flexible phenomenological model ( |l6|) , based on 
two exponential decays, appears more promising. In this 
case the best fltted lineshapes (quantum corrected and 
resolution convoluted), including the thermal relaxation 
(using the Q — > values of specific heat ratio and ther- 
mal diffusivity from [^3|) through Eq. (p^, are again 
reported in Figs. Q and |5|, and a marked improvement of 
the agreement is clearly visible. 

In Fig. ^ we also report the classical and resolution de- 
convoluted lineshapes of S{Q, uj) and Cl {Q, ^jj) as builded 
from the fitting parameters, in the region Q < 14 nm~^, 
the most significant as far as the positive dispersion task 
is concerned. The maxima of the latter function have 
been utilized to determine the apparent sound velocity 
of the system. 

Let us now discuss the different parameters entering 
Eq. (p^), starting from the relaxation times: in Fig. 
^ we report the wavevector dependence of such quanti- 
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ties at both temperatures. The presence of the double 
timescale is confirmed by the difference by about one or- 
der of magnitude between the two parameters Tq, and . 
Concerning Tq it is worth to point out the non negUgible 
effect of resolution: the value of 1/ HWHM is compara- 
ble to the relaxation time itself so that it is expected to 
affect this fitting parameter, particularly at low Q. As 
a consequence it is unfortunately not possible to inspect 
any temperature dependence of Ta (Q,T). Completely dif- 
ferent is the situation as far as the fast timescale is con- 
cerned. The determination of r^^ is expected to be defi- 
nitely more reliable and on the basis of our analysis the 
fast dynamics turns out to be temperature independent 
at the investigated Q-values. The influence of these re- 
laxation mechanisms on the dynamics can be inferred 
looking at the uj{Q)t{Q) values as done before for the 
thermal channel: as far as the slow mechanism is con- 
cerned the system response turns out to be "solid like" 
across the full first Brillouin pseudo-zone (see the insets 
in Fig. 0a). The fast relaxation accounts instead for a 
moderately viscous response: uj{Q)t{Q) ~ 1 around the 
peak value at Q ~ 8 nm~^and remains quite lower all 
around (insets in Fig. ^). Following the arguments of 
Sec. |II B| , we expect the slow relaxation to be responsi- 
ble for the narrow quasi elastic contribution and for the 
position of the inelastic features of S{Q, oj) while the fast 
process should drive mainly the acoustic damping (Bril- 
louin linewidth) also accounting for the broader part of 
the Mountain peak |Q. 

To go further and point out the real quantitative in- 
fluence of such two channels is necessary to look at the 
strength parameters, or better at the rate a{Q). Such 
quantity, defined as the ratio A^/ (A^ -I- A^) , has its 
maximum of ~ 0.20 at low Q and decreases at increasing 
the wavevector as shown in Fig. ^. This means that the 
fastest process plays the major role in characterizing the 
dynamics of liquid lithium all over the investigated first 
Brillouin pseudo-zone. Such statement can be better un- 
derstood looking at the apparent sound speed behavior 
as plotted in Fig. ||. The sound velocity moves from the 
low Q value (isothermal in the IXS window) to its high 
frequency value Coo(Q). As previously noted, such behav- 
ior has been observed in several MD simulations of alkali 
metals and in some cases experimentally, even if with 
some ambiguity due to kinematic restrictions/incoherent 
contribution on neutron side and resolution broadening 
in the pioneer IXS experiments . What is less clear is 
the microscopic origin of such positive dispersion. On the 
basis of our analysis, once unambiguously found out an 
experimental evidence of a double timescale, we are able 
to ascribe the positive dispersion of the sound velocity to 
the faster process. We reported in fact in Fig. |9|the value 
of the solid- like response associated to the single slow pro- 
cess and to both slow and fast relaxation processes. As 
expected the "jump" in sound velocity associated to the 



slow a-process ( ^Wq((5) + A^((5) — ujo{Q) /Q) has al- 
ready occurred at the lower IXS Q values, but it is quanti- 
tatively negligible as compared to the full positive disper- 
sion of the system. Despite of the fact that the condition 
'-^iiQ)Tti{Q) S> 1 is never satisfied, the faster process, due 
to its dominant strength, accounts almost entirely for the 
viscous- elastic transition. As can be seen in Fig. ^, the 
apparent sound velocity Capp — uJi{Q)/Q, never reaches 
the infinite frequency sound velocity deduced by the fit, 

coo(g) = y/c.2(g) + A2(g) + A2(Q)/g. 

At first sight this discrepancy may be explained on the 
fact that uJi{Q)Tf^{Q) is always 1. However Coo{Q) can 
be also evaluated theoretically starting from the knowl- 
edge of the interatomic potential and the pair distribu- 
tion function (see for example ref. 0). In Fig. ^we also 
report the values of the unrelaxed sound speed c^((9) 
computed in this way, they turns out to be lower than 
those determined from the fitting of the spectra at both 
the temperatures. Moreover, c*^{Q) is actually reached 
by the apparent sound speed at Q = 9 nm^^, i.e. when 

A possible explanation of the inconsistency between 
c^(g) and Coo(Q) - ^u;^,{Q) + Al{Q) + Al{Q)/Q can 
be related to the arbitrary choice of the memory function 
details. Even though the presence of a two timescale 
dynamics is directly seen in the raw experimental data 
(particularly in the 7 < Q < 14 nm~^ region, the quasi 
elastic peak reflects the superposition of two peaks of 
different widths), the detailed shape of Ml{Q, t) can have 
in principle more complicated features than the double 
exponential ansatz of Eq. (|l6|). In particular we already 
mentioned the major drawback of such an assumption, 
namely the cusp at i = 0. It is reasonable to think that an 
exponential decay forced to represent a more complicated 
time dependence can give an accurate estimate as far as 
the r is concerned, while at short times the lack of a zero 
second derivative inescapably leads to an overestimate of 
Ml{Q, t = 0), i.e. just the positive dispersion amplitude: 
a similar effect, i.e. the overestimation of the Coo{Q 
0) deduced by the fit of BLS data using instantaneous 
approximation for the microscopic process with respect 
to the experimental IXS current spectra maxima, has 
been recently shown in polybutadiene . 

The last important aspect concerns the relative weights 
of the two processes and how to relate them to the acous- 
tic attenuation. The role of the thermal effects has been 
already discussed. Looking at the generalized viscosity 
contributions A^r^ and A^Tq (see Fig. |l^) we find that 
they are nearly comparable, (the strength differences are 
nearly balanced by the different r), but once again only 
the microscopic (faster) relaxation process can be related 
to the Brillouin linewidth being the only one close to the 
ujT << 1 condition in the explored g-range. Moreover, in 
the spirit of the hydrodynamic generalization, the sum of 
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such contributions is expected to have a dependence. 
Such prediction is quite well reproduced by our fits: only 
at very low Q a slight deviation occurs, ascribable to the 
slow part; once again the reason may be an excess esti- 
mate of the Of— relaxation time at small Q. Finally, in 
Fig. ^ the longitudinal viscosity as outcome of the fit 
(the full area under M{Q,t)) is compared to an exper- 
imental determination of shear static {Q — > 0) viscosity 
p3[ . Noting that in molten alkali metals the shear part 
usually accounts for a 30 — 40% of the longitudinal part, 
the agreement is satisfactory, except for the very low Q 
points where finite resolution function effects prevent an 
accurate check. 



be reached in the very low Q (adiabatic) and high Q (free 
particle) 

After an initial nearly linear dispersion, structural 
effects take place suppressing the sound propagation 
around Qm/'^ due to strong negative interference. With 
increasing Q values, the points in Fig. |I2| show not only 
a second pseudo-BZ, but also a series of oscillations that 
damp out with increasing Q - here, loi{Q) is approach- 
ing the single particle behavior. These oscillations are 
in anti-phase with those of S{Q) (see Fig. ^) and can 
therefore be associated with the local order in the liquid. 

V. CONCLUSION 



B. The "high Q region" 

Above Qrm collective excitations are still evident in 
I{Q,Lu) for nearly one more Brillouin pseudo-zone, even 
if the analysis of the previous section is not wholly sat- 
isfactory because of the worsening of the experimental 
resolution. Moreover, single particle features gradually 
emerge giving rise to a crossover between two completely 
different dynamical regimes. As the density fluctuation 
wavelength matches the mean interparticle distance, a 
transition between strongly correlated to incoherent mo- 
tion occurs. The lineshape associated to a quantum free 
particle is well known to be Gaussian. Once that Sq{Q) 
is determined as before from the first two moments sum 
rules, there are no unknown parameters. While in the 
previous section our sample has been supposed to be 
single-component, we must now, due to the emergency 
of additional quantum effect, be more accurate by taking 
into account isotopic effects. From the natural abun- 
dance ratios (92% of ^Li, 8% of ^Li) one finds that 

5,(Q,..) = -^^^e-(--)V2-? (22) 

2—6,7 

where Ci are the above mentioned isotope concentrations 
while LUi — hQ^/2mi and af —KbTQ^ /rrii are the first 
(central) and the second (relative to Ui) frequency mo- 
ments of the two isotopes. In Fig. || the experimental 
spectra are reported together with the theoretical line- 
shapes (resolution convoluted) of Eq. ( ^ ) in order to 
emphasize all the crossover features. To be more quan- 
titative we report in Fig. |l2| the Stokes current spectra 
maxima as function of the exchanged wavevector. In- 
deed, in this "high Q" region (above (5m=24 nm~^) this 
quantity has been calculated starting from the spectra of 
Fig. H (we expect the effect of the resolution broadening 
to be negligible at such high Q' s). For this reason, due 
to the quantum features of the S{Q,uj), the quantity oji 
takes different values in the Stokes/ Antistokes sides of 
the spectra Q. The asymptotic values are expected to 



In the present work inelastic X-ray scattering has been 
utilized to study in detail the main features of the col- 
lective dynamics in liquid lithium, a system where the 
use of inelastic neutron scattering is quite difficult and 
basically impossible in the low Q region. This IXS ex- 
periments have allowed for a precise assessment of the 
different relaxation processes which affect the spectral 
shape of the dynamic structure factor. Specifically, the 
quality of the data required an analysis where one has to 
invoke the simultaneous presence of three distinct decay 
channels of the collective memory function. The first one 
is associated with the coupling between density and tem- 
perature fluctuations. As in other alkali metals, although 
its relevance is rather small, this process cannot be ne- 
glected and plays a different role in respect to ordinary 
non-conducting liquids. Much more important is the un- 
ambiguous evidence of two well separated timescales in 
the decay of the ML{Q,t) - the portion of the memory 
function directly associated with generalized viscosity ef- 
fects. In this respect, our IXS data prove that the tradi- 
tional use of a single time-scale (viscoelastic model) can- 
not account for the detailed shape of S{Q,uj), and that 
it only provides a guidance for a qualitative description 
of different dynamical regimes. 

Despite the effect of an additional time-scale to the vis- 
coelastic model has been tested on MD results in the case 
of Rubidium (see pag. 248 of ref. 0), and theoretical in- 
terpretations of such framework (in particular devoted 
to frame the slow process within Mode Coupling theory) 
have been recently given p3| we give in this work, to the 
best of our knowledge, the first experimental evidence of 
the necessity of such an assumption to correctly repro- 
duce the measured S{Q, oj), stressing the role of the con- 
tributions from the different relaxations to the S{Q,lu) 
in the considered Q — E regions. The presence of a two 
time-scales mechanism leaves open, however, a point of 
crucial importance: namely, their physical origin and, 
most important, of the fast one. 

The slow process, referring to the terminology used 
to describe glass forming systems, is related to the a- 
relaxation responsible for the liquid-glass transition in 
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systems capable to sustain strong supercooling. The 
same process can be framed within kinetic theory in 
terms of " correlated collisions" : in mode coupling ap- 
proaches, the onset of these correlation effects is traced 
back to the coupling to slowly relaxing dynamical vari- 
ables, and specifically in the liquid region, to long lasting 
density fluctuations. The corresponding decay mecha- 
nism (often referred to as "structural relaxation"), al- 
most negligible in the initial coUisional region, controls 
the dynamics at intermediate and long times. For Q — > 
the amplitude of the slow portion of the memory func- 
tion turns out to vanish as Q^, whereas at increasing 
wavevectors it becomes proportional to Q . At suffi- 
ciently large Q, the relevance of these couplings as long 
lasting decay channel decreases markedly. In fact, in a 
rather wide portion of the explored range of wavevectors, 
we find that the amplitude /\\{Q) of the total memory 
function behaves approximately as Q^. Therefore, outside 
the hydrodynamic region, the normalized weight of the 
slowly varying contributions to ML(Q,t) approximately 
decreases proportionally to Q/Q^ = l/Q- This weight 
may reasonably be identified with the dimensionless am- 



plitude factor a{Q) in the model of Eq. (16). A quantita- 
tive understanding of the " slow" timescale would require 
a full evaluation of the mode-coupling contribution at 
different wavevectors, as in fact done in the above men- 
tioned MD study |||. 

The situation is much more involved as far as the fast 
process is concerned and, in fact, its microscopic inter- 
pretation is still matter of debate. In the present work 
we have given an experimental evidence that, in liquid Li, 
the fast process is indeed the dominant one, and it mainly 
controls both the sound speed dispersion and attenuation 
i.e. the position and the width of the Brillouin compo- 
nent of the spectra. This observation justifies the need 
for a detailed discussion on the nature of the fast relax- 
ation channel. Within generalized kinetic theory the fast 
initial decay of both single particle and collective mem- 
ory function is traced back to coUisional events, which 
are short ranged both in space and in time. In partic- 
ular, in the collective case, the short timescale l/7i(Q) 
can be associated with the duration of a rapid structural 
rearrangement occurring over a spatial range = 2'K/Qm- 
In most liquids, at low and intermediate wavevectors, the 
typical value of the time l/ji{Q) is of the order of a few 
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s, with some structural oscillations of decreasing 



trend |3S|]. At Q > Qm, the structural effects damp out, 
the decrease becomes more marked, and, when the single- 
particle aspects start to prevail, l/ji{Q) oc Q^^. An in- 
spection of the results reported in a recent MD study in 
liquid lithium near melting basically confirms these 
qualitative features as well as the actual values of the 
"short time" found in our data analysis (see Fig. |^). 
Although the description of the fast process in term of 
interparticle collision is a possible way to account for the 
dynamical features at short times, it does not allow a 



deep insight into the physics behind it. 

Although the description of the fast process in terms of 
local coUisional events has been widely used in the past to 
account for the short-time dynamics, other mechanisms 
cannot be excluded. We have already mentioned in sec. 
II that a formally similar "microscopic mechanism" has 
been introduced in the analysis of the Brillouin spectra 
of glassformers |25| ; more recently, an analogous micro- 
scopic relaxation process has been detected in a simulated 
harmonic glass |^^. Generally speaking, the inclusion of 
a relaxation process in a memory function reflects the 
existence of decay channels by which the energy stored 
in a given degree of freedom ("mode") relaxes toward 
other modes. The problem is then the physical identifl- 
cation of the modes entering the fast relaxation process, 
and the ultimate origin of the latter in a widely different 
class of systems ranging from fluids to glasses. In gen- 
eralized kinetic theories the mechanism is traced back to 
a mainfold of non- hydrodynamic, phase-space "kinetic" 
modes, whose non-conserved character ensures a rapid 
decay of the memory function |39|j4l| ]. A different ap- 
proach to describe the fast part of the memory function 
in disordered systems relies on the normal modes ("in- 
stantaneous" in normal liquids) analysis of the atomic 
dynamics . In this case one uses a framework (dynam- 
ical matrix, etc.) formally similar to the one customarily 
adopted for harmonic crystals; however, owing to the lack 
of translational symmetry of the system, it turns out ]4C| ] 
that the eigenstates cannot anymore be represented by 
planewaves (PW) even at relatively small wavevectors 
p3[ . As a result, when probing the system at a spe- 
cific wavevector Q one effectively detects a transfer from 
a PW mode toward other PWs of different wavevectors. 
Assuming that it is actually this energy flow which causes 
the fast relaxation, we deal with a mechanism whose ul- 
timate origin is the topological disorder, and not a truly 
dynamical event. Ordinary liquids (such as the one con- 
sidered here) are certainly "disordered systems" and at 
the high frequency considered here (cjTq >> 1), they 
can be considered as "frozen"; thereforethe description 
used for glasses can be applied as well. At the present 
stage of the theory our experimental data cannot sup- 
port either interpretation, or not even ascertain whether 
the distinction is largely semanthic, or not. However, the 
posssibility that the phenomenology reported here for a 
simple liquid can be understood within the same men- 
tal framework developed for more complex liquids and 
glasses, is certainly interesting. Further studies on other 
simple fluids, as for example neon seem to suggest 
that is indeed the case. 
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"Density Fluctuations in Molten Lithium..." by T.Scopigno et al. 




Frequency shift (ps ^) 

FIG. 1. IXS spectra at selected Q < Qm values, in the (9 9 9) configuration, at two different temperatures. The full line is 
the best fit taking into account thermal, slow and fast relaxations. 
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"Density Fluctuations in Molten Lithium..." by T.Scopigno et al. 




FIG. 2. Comparison between the S{Q) experimentally determined as discussed in the text and trough MD simulations ( ) 

of ref.[28]. Triangle (A): (11 11 11) data, Dots: (o) (9 9 9) data, Lozenges {<}): (7 7 7) data. The fuh symbols are the values 
corrected for asymmetry resolution effects, at the higher Q's such correction is negligible and is not reported. The inset shows 
on an enlarged scale the low Q region. 
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FIG. 3. IXS spectra at Q values above the main peak of the structure factor, in the (7 7 7) configuration, at the lower 
temperature. The full line is the theoretical lineshape expected in the single particle, high Q, limit. 
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"Density Fluctuations in Molten Lithium..." by T.Scopigno et al. 




FIG. 4. Comparison between the value of loq{Q) as calculated from its definition ( ) and from the three times fit (A) 

leaving it as a free parameter. The normalized second frequency moment of l{Q,u}) is also reported (•) to emphasize resolution 
and quantum effects. In the inset we show the same non-normalized quantities. 
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"Density Fluctuations in Molten Lithium..." by T.Scopigno et al. 




FIG. 5. Inelastic x-ray scattering spectrum of Lithium at T = 475 K and Q = 7.1 nm~^ Dots (•): experimental data; Dashed 
line ( ): simple hydrodynamic fit (DHO); Dotted lines ( ): simple viscoelastic model (thermal and one relaxation 

time) showing two different sessions of comparable x'^j Full line ( ): thermal and two relaxation times, structural and 

microscopic. The inset shows the mean square deviation of the different models. 



18 



"Density Fluctuations in Molten Lithium..." by T.Scopigno et al. 



c/3 

o 



3 




Q = 1.4 mn 2.0 



40 60 



s 



o 

C/3 



3 

3 



CO (ps" ) 



FIG. 6. Dynamical structure factor S{q,u}) 



-) and longitudinal current correlation function Cl{Q,u}) ( ) in the 



low-Q region as obtained from the two times fitting model (classical and resolution deconvoluted) . The (symmetical) peak 
positions of the latter quantity have been utilized to determine the apparent sound speed of the system. 



19 



"Density Fluctuations in Molten Lithium..." by T.Scopigno et al. 



C/3 



O 



0.06 



0.04 



■5: 



1.0 - 



0.5 



0.0 



5 . -o 



10 



20 



0.02 - 



• • s „ Q 

• Q 



O 0.00 



a 



20 



15 



I 10 



0.5 - 



V 



08 



• o 



o8o«» • ' • •• 



. _i le 

Q(nm-') 



0.0 







10 15 20 25 



Q (nm ) 



FIG. 7. Relaxation times at different temperatures from the fits. Full symbol: T = 475 K, open symbol T = 600 K. The 

timescale corresponding to the experimental resolution (3.0 meV) is also reported ( ) to show how it could affect the lower 

relaxation time determination. 
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FIG. 8. Ratio between the slow and the total relaxation strength. Full dots (•) T = 475 K, open dots (o) T = 600 K. 
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FIG. 9. Apparent sound speed (o), as from the maxima of J(Q,uj) build from the fit parameters (therefore resolution 

deconvoluted) , compared to the hydrodynamic limit of ref. [33] (<■ — arrow) and to the high frequency limit c^(Q) ( ) of 

ref. [28]. The following quantities, deduced from fitting parameters, are also reported: the isothermal speed of sound co(<5) 
(0), the high frequency limit associated to the slow process Vwg(Q) + A2(Q)/Q (-V-) and finally the total Coc(Q) (- A -) 
with all the strenghts included. The full lines/symbols denote T = 475 K data, the dotted lines/full symbols are relative to 
T = 600 K. 
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FIG. 11. Longitudinal viscosity from the fit (A„ra+A^r^)/(3^ compared to shear viscosity data of ref. The full lines/symbols 
denote T = 475 K data, the dotted lines/open symbols are relative to T = 600 K. 
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FIG. 12. Maxima of the current spectra (•) in the full range of explored wavevectors at T = 475 K, compared to the 

asiritotic opposite behaviors in the low Q (hydrodynamic ( ) of ref. [33]) and high Q (quantum free particle ( ), see 

text) limits. As far as (7 7 7) data are concerned we reported the rough current maxima as calculated from the spectra. To 
estimate resolution effect we also include in this case the high Q limit value as deduced by the free particle model folded with 
instrumental resolution. 
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